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Abstract. Motivated by juggling sequences and bubble sort, we examine per- 
mutations on the set {1, 2, . . . , n} with d descents and maximum drop size k. 
We give explicit formulas for enumerating such permutations for given integers 
, k and d. We also derive the related generating functions and prove unimodality 

1 and symmetry of the coefficients. 
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1. Introduction 

There have been extensive studies of various statistics on S n , the set of all 
permutations of {1,2, ... ,n}. For a permutation tt in S n , we say that tt has a 



> 



drop at i if 7Ti < i and that the drop size is i — 7T;. We say that tt has a descent at 
i if TTi > One of the earliest results [5] in permutation statistics states that 

i^h ' the number of permutations in S n with k drops equals the number of permutations 

with k descents. A concept closely related to drops is that of excedances, which is 
just a drop of the inverse permutation. In this paper we focus on drops instead of 
excedances because of their connection with our motivating applications concerning 
bubble sort and juggling sequences. 

Other statistics on a permutation tt include such things as the number of in- 
' versions, that is, |{(i,i) : i < j, TTi > and the major index of tt (i.e., the 

I/"") , sum of i for which a descent occurs) . The enumeration of and generating functions 

for these statistics can be traced back to the work of Rodrigues in 1839 [TU] but 
was mainly influenced by McMahon's treatise in 1915 [5]. There is an extensive 
literature studying the distribution of the above statistics and their g-analogs (see 
for example Foata and Han [4] or the papers of Shareshian and Wachs [Til E2] for 
more recent developments). As noted above, the drop statistic that we study is 
closely related to the excedances statistic. The distribution of the bivariate statis- 
tics (descents, excedances) can be found in Foata and Han [51 equations (1.15) and 
(1.16)]. _ 

' This joint work originated from its connection with a paper [2] on sequences that 

can be translated into juggling patterns. The set of juggling sequences of period 
n containing a specific state, called the ground state, corresponds to the set B n ,k 
of permutations in S n with drops of size at most k. As it turns out, B n ^ can also 
be associated with the set of permutations that can be sorted by k operations of 
bubble sort. These connections will be further described in the next section. We 
note that the maxdrop statistic has not been treated in the literature as extensively 
as many other statistics in permutations. As far as we know, this is the first time 
that the distribution of descents with respect to maxdrop has been determined. 

First we give some definitions concerning the statistics and polynomials that 
we examine. Given a permutation tt in S n , let Des(7r) denote the descent set, 
{1 < i < n : TTi > 7Ti+i}) °f tt and let des(7r) = | Des(7r)| be the number of descents. 
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B n ,k{x) ='^b n , k {r)x r = E ;c- 



We use maxdrop(7r) to denote the value of the maximum drop (or maxdrop) of tt, 

maxdrop(7r) = max{ i — n(i) : 1 < i < n }. 

Let B n .k — {tt G S n : maxdrop(7r) < k}. It is known, and also easy to show, that 
\B n ,k\ = fc!(fc + e.g., see Thm. 1] or 8, p. 108]. Let 

b n ,k(r) = |{> e B n<k : des(?r) = r}\, 

and define the (fc-maxdrop-restricted) descent polynomial 

r des(7i 

■r>0 JreS n ,n, 

Examining the case of k = 2, we discovered that the coefficients b n ^{r) of B n ^{x) 
appear to be given by every third coefficient of the simple polynomial 

(l + x a )(l + x + x 2 ) n - 1 . 

Looking at the next two cases, k = 3 and k = 4, yielded more mysterious polyno- 
mials: b n> s(r) appeared to be every fourth coefficient of 

(1 + x 2 + 2x 3 + x A + x 6 )(l + x + x 2 + x 3 ) n - 2 

and 6,1,4 (r) every fifth coefficient of 

(1 + x 2 + 2x 3 + Ax A + 4x 5 + 4a; 7 + 4x 8 + 2x 9 + x w + x 12 )(l + x + x 2 + x 3 + a? 4 )™" 3 . 

After a fierce battle with these polynomials, we were able to show that & n fc(r) is 
the coefficient of 

u r(k+l) m 

the polynomial 

P k (u)(l + u + ---+u k ) n ' k (1) 

where 

k fe y , x 

fl b («)=x;^-i(«* +i )(«* +i -i) i E (*)«-*, (2) 

and Afc denotes the fcth Eulerian polynomial (defined in the next section). Further 
to this, we give an expression for the generating function Bfc(z, y) = Yl n >o B n ,k(y)z n , 
namely 



i + e (mv) - E ( h t 1 ) (y - w-^t-iiv) ) z 



B k (z,y) = — ' 1 



k+1 'Jfe + 1 



i=i 



We also give some alternative formulations for which lead to some identities 
involving Eulerian numbers as well as proving the symmetry and unimodality of 
the polynomials B n ^{x). 

Many questions remain. For example, is there a more natural bijective proof for 
the formulas that we have derived for B nk and B&? Why do permutations that 
are fc-bubble sortablc define the aforementioned juggling sequences? 



2. Descent polynomials, bubble sort and juggling sequences 
We first state some standard notation. The polynomial 

^,des(Tr) 



A n (x) = E x 



is called the nth Eulerian polynomial. For instance, Ao(x) — A\(x) — 1 and 
A^ix) = 1 + x. Note that B n ^{x) = A n (x) for fc > n— 1, since maxdrop(7r) < n— 1 
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for all 7r € S n . The coefficient of x k in A n (x) is denoted (?) and is called an 
Eulerian number. It is well known that ([B]) 

1 — w \ / n\ h z n 



,(«;-l)2_ lt; \fc 
fc,n>0 



The Eulerian numbers are also known to be given explicitly as ([3JE]) 

We define the operator bubble which acts recursively on permutations via 

bubb\e{LnR) = bubble(L)i?n. 

In other words, to apply bubble to a permutation ir in <S„, we split tt into (possibly 
empty) blocks L and R to the left and right, respectively, of the largest element 
of 7r (which initially is n), interchange n and R, and then recursively apply this 
procedure to L. We will use the convention that bubble(0) = 0; here denotes the 
empty permutation. This operator corresponds to one pass of the classical bubble 
sort operation. Several interesting results on the analysis of bubble sort can be 
found in Knuth 8, pp. 106-110]. We define the bubble sort complexity of 7r as 

bsc(7r) = minjfc : bubble fe (7r) = id}, 

the number of times bubble must be applied to ir to give the identity permutation. 
The following lemma is easy to prove using induction. 

Lemma 1. (i) For all permutations ir we have maxdrop(7r) = bsc(7r). 
(ii) The bubble sort operator maps B n _k to B n _k-i- 

The analysis of algorithms similar to bubble sort has been instrumental in gen- 
erating interesting research. For example, the analysis of stack sort in Knuth.0 pp. 
242-243] gave rise to the area of pattern avoiding permutations. The stack sort 
operator stack is defined by stack(Lni?) = stack(i) stack(i?)n. We see below that 
stack sort is at least as efficient as bubble sort. 

Lemma 2. For all ir G S n , if bubble fe (7r) = id then stack fc (7r) = id. 

Proof. The proof of Lemma [5] follows from the following claim: 

If A = a\02 ■ ■ ■ a n = LmR is a sequence of distinct positive integers and m = 
maxiOi, then either maxdrop(j4) = l — a\ or maxdrop(stack(A)) < maxdrop(A) — 1. 

The Claim is certainly true for n = 1. Suppose the claim is true for n' < n. If 
maxdrop(A) = 1 — oi, we are done. We may assume that maxdrop(A) > 1 — a%. 
This implies that the maxdrop of A does not occur at the entry where m is located. 
For the ith entry in stack(i), the maxdrop of stack(L) at i is reduced by one by 
induction. For the jth entry in R, the maxdrop of the corresponding entry in 
stack(A) is reduced by 1. Thus, the claim is proved by induction. □ 

The class of permutations B n ^ appears in a recent paper [2] on enumerating 
juggling patterns that are usually called siteswaps by (mathematically inclined) 
jugglers. Suppose a juggler throws a ball at time i so that the ball will be in the air 
for a time t{ before landing at time U + i. Instead of an infinite sequence, we will 
consider periodic patterns, denoted by T = (ti,t%, . . . , t n ). A juggling sequence is 
just one in which two balls never land at the same time. It is not hard to show [1] 
that a necessary and sufficient condition for a sequence to be a juggling sequence 
is that all the values U+i (mod n) are distinct. In particular, it follows that that 
the average of ti is just the numbers of balls being juggled. Here is an example: 
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If T = (3, 5, 0, 2, 0) then at time 1 a ball is thrown that will land at time 1 + 3 = 4. 
At time 2 a ball is thrown that will land at time 2 + 5 = 7. At time 3 a ball is 
thrown that will land at time 3 + = 3. Alternatively one can say that no ball is 
thrown at time 3. This is represented in the following diagram. 




1 2 3 4 5 6 7 



Repeating this for all intervals of length 5 gives 











1 2 3 4 5 


6 7 8 9 10 



For a given juggling sequence, it is often possible to further decompose into 
shorter juggling sequences, called primitive juggling sequences, which themselves 
cannot be further decomposed. These primitive juggling sequences act as basic 
building blocks for juggling sequences [2J. However, in the other direction, it is 
not always possible to combine primitive juggling sequences into a longer juggling 
sequence. Nevertheless, if primitive juggling sequences share a common state (which 
one can think of as a landing schedule), we then can combine them to form a 
longer and more complicated juggling sequences. In [2J, primitive juggling sequences 
associated with a specified state are enumerated. Here we mention the related fact 
concerning B n ^\ 

There is a bijection mapping permutations in B n ,k to primitive juggling sequences 
of period n with k balls that all share a certain state, called the ground state. 

The bijection maps 7r to (f>(ir) = (t±, . . . , t n ) with = k — i+Tti. As a consequence 
of the above fact and Lemma [TJ we can use bubble sort to transform a juggling 
sequence using k balls to a juggling sequence using k — 1 balls. 

To make this more precise, let T = {t\, . . . ,t n ) be a juggling sequence that 
corresponds to 7r 6 B n ,k, and suppose that T" = (si, . . . , s n ) is the juggling sequence 
that corresponds to bubble(7r). Assume that the ball B thrown at time j is the one 
that lands latest out of all the n throws. In other words, tj +j is the largest element 
in {ti+i}2 =1 - Now, write T — LtjR where L — (t\, . . . , tj-i) and R = (tj+i ,...,£„). 
Then we have 

T' = f k (T) = f k (L)Rs, 

where s = tj +j — (n + 1). In other words, we have removed the ball B thrown at 
time j and thus throw all balls after time j one time unit sooner. Then at time n we 
throw the ball B so that it lands one time unit sooner than it would have originally 
landed. Then we repeat this procedure to all the balls thrown before time j. 

3. The polynomials B n ^ k {y) 

In this section we will characterise the polynomials B n j;(y). This is done by 
first finding a recurrence for the polynomials and then solving the recurrence by 
exploiting some aspects of their associated characteristic polynomials. The latter 
step is quite involved and so we present the special case dealing with B n ^{y) first. 

3.1. Deriving the recurrence for B n k- We will derive the following recurrence 
for B n , k (y). 
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Theorem 1. For n > 0, 



B„+k+l,k 



M-gCt 1 )"- 1 ' 



i-1 



B n +k+l-i,k{y) 



(4) 



with the initial conditions 



Bi, k (y)=Mv)> 0<i<k. 



We use the notation [a, 6] = {i e Z : a < i < 6} and [b] = [1,6]. Let A = 
{ai, . . . , a„} with a\ < ■ ■ ■ < a n be any finite subset of N. The standardization of a 
permutation ir on A is the permutation st(7r) on [n] obtained from 7r by replacing the 
integer a* with the integer i. Thus it and st(7r) are order isomorphic. For example, 
st(19452) = 15342. If the set A is fixed, the inverse of the standardization map is 
well defined, and we denote it by st^ (<r); for instance, with A = {1,2,4,5,9}, we 
have stT 1 (15342) = 19452. Note that st and stT 1 each preserve the descent set. 

For any set S C [n — 1] we define A n ,k(S) — {n e B n ,k ■ Dcs(7r) D S} and 



Note that t n (S) = in the case that n — 1 is not a member of 5*. Now, for any 
permutation 7r = 7Ti . . . 7r„ in A n ,k(S) define 

f(n) = (o;X), where a = st(ni . ..n n -i-i),X = {7r„_;, . . . ,n n } and i = t n (S). 

Example 1. Let S — {3,7,8}, and choose the permutation tt = 138425976 in 
A 9t3 (S). Notice that Dcs(tt) = {3,4,7,8} D S. Now t g (S) = 2. This gives /(tt) = 
(a,X) where a = st(138425) = 136425 and X = {tt 7 , tt 8 , m,} = {6,7,9}. Hence 
/(138425976) = (136425, {6, 7, 9}). 

Lemma 3. For any tt in A n ,k(S), the image f(ir) is in the Cartesian product 



Proof. Given tt € A n ,k(S), let /(7r) = (<x, X). Suppose i = t n (S). Then there are 
descents at positions n — i,...,n—l (this is an empty sequence in case i = 0). Thus 



where the last inequality follows from the assumption that maxdrop(7r) < k. Hence 
X is an (i + l)-element subset of [n — k, n], as claimed. Clearly a e S n -i-i- 

Next we shall show that o is in A n -i-i.k- Notice that the entries of (iri, . . . , 7r„_i_i) 
that do not change under standardization are those irg which are < 7r„ . Since these 
values remain unchanged, the values t — ire are also unchanged and are thus < k. 

Let (7r a (i): • • • j ^afm)) be the subsequence of values which are > 7r„. The smallest 
value that any of these may take after standardization is n n > n — k. So c a (j) > 
7r„ > n— k for all j e [1, m]. Thus a(j) — a a ^ < a(j) — (n—k) = k— (n — a(j)) < k 
for all j £ [1, m]. Therefore £ — og < k for all I € [l,n — i — 1] and so o e A n -i-\,k- 

The descent set is preserved under standardization, and consequently a is in 
A n -i-i,k(S n [n — i — 2]), as claimed. □ 

We now define a function g which will be shown to be the inverse of /. Let ir 
be a permutation in A m ,k(T), where T is a subset of [m — 1]. We will add i + 1 
elements to tt to yield a new permutation a in A m +i.k(T U [m + 1, m + Choose 
any (i + l)-element subset X of the interval [m + i + 1 — k, m + i + 1], and let us 
write X = {x\, . . . ,Xi + i}, where x\ < ■ ■ ■ < x i+1 . Define 

g(n, X) = sty 1 (7Ti . . . 7r m ) x i+1 Xi ...xi, where V = [m + i + 1] \ X. 



t n (S) = max{i eN:[n-!,n-l]CS}. 




where (^) denotes that set of all m-element subsets of the set X. 



n > lT n -i > 7T. 



n-i+l > ■ ■ ■ > 7T„_i > 7T„ > n - k, 
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Example 2. LetT = {1}, and choose the permutation ir = 3142 inAi^T). Notice 
that Des(7r) = {1, 3} D T. Choose i = 2 and select a subset X from [4 + 2 + 1 — 
3, 4 + 2 + 1] = {4, 5, 6, 7} of size i + 1 = 3. Let us select X = {4, 6, 7}. Now we 
have g(n,X) = sty 1 (3142) 764 = 3152764, where V is the set [4 + 2 + 1] \ {4, 6, 7} = 
{1,2,3,5}. 

Lemma 4. If (ir, X) is in the Cartesian product 

'[m + i + l — k, m + i + 1] 



A m , k ( T ) x . 

i + 1 

for some i > then g(ir, X) is in 

A m +i+i,k(TU [m + l,m + i}). 

Proof. Let a = g{ir,X). For the first m elements of a, since <jj > irj for all 
1 < j < m, we have j — aj < j — nj which gives 

maxjj — Uj : j G [m]} < maxjj — TTj : j & [m]} < k. 

The final i + 1 elements of <r are decreasing so the maxdrop of these elements will 
be the maxdrop of the final element, 

m + i + 1 — (Tm+i+i = m + i + \ — x\ < m + i + 1 — (m + i + 1 — k) = k. 

Thus maxdrop(er) < k and so a £ B m +i+i, k . The descents of a will be in the set 
TU[m+l, m+i] since descents arc preserved under standardization and the final i+1 
elements of a are listed in decreasing order. Hence o e A m+i+ i^(TU [m + l,m+i]), 
as claimed. □ 

Lemma 5. The function f is a bijection, and g is its inverse. 

Proof. Given any (a,X) e x A m ,k( T ) whcrc T <Z [m - 1], let 

it = g(cr,X). We have 

n = st [m 1 +j+i]\x( cr i • ■■<r m )xj + iXj ...xi e An+j+i,fe(ru [m+ l,m + j]), 

where X = {xi, . . . and x\ < ■ ■ ■ < Xj + \. Let S = Tu[m+l,m+j]. Clearly 

Des(g{a,X)) D S and i = t m+j+1 {S) = j. So f(g(a,X)) = (t,Y) whcrc 

Y = {n m+ i, . . . ,7r m+ j+i} = {xi, . . . , Xj + i} = X 

and 

r = st(st^ ro+ _. +1]VX ((7i . . . a m )) = a. 

Hence f(g(a,X)) = {a,X). Given tt e A n . k {S), let /(tt) = (cr,X) with X = 
{xi, . . . , Xi + i} and cr = st(7Ti . . . 7r n _(j + i)). We have 

3((T, X) = St^ n]xx ((Ti . . . (Tn-i-^Xi+x . . . X\ 

= ^%n\\X (stC 71 "! ' ' • 7r n-i-l)) 7r n-J ■ • ■ 7T„ 
= 7Tl . . . TTn— i— 1 ^n — 2 ■ ■ • 7T n 
= 7T. 

Hence g(f(n)) = n. □ 
Corollary 1. Le£ a n _k{S) = \A n ,k{S)\ and i = t n (S). Then 

(k + 1\ 

an,fc(5) = ( . + 1 Ja„_(i + i), fe (5n[l,n- (i + 1)]). 
Proposition 1. For a/Z n > 0, 

BnAv + 1) = E ( fc 1 + 



DESCENT POLYNOMIALS FOR PERMUTATIONS WITH BOUNDED DROP 



7 



Proof. Notice that 

B n , k (y + 1)= J2 (2/ + l) d ° sW 



E E 



= E E y |s| 

7reB n , fc SCDes(Tr) 

= E ^' 5| E 1= E 2/ |s| «^(^)- 

SC[n-l] 7rG.A„,fc(S) SC[n~l] 

From Corollary [TJ multiply both sides by y\ s \ and sum over all S C [n — 1]. We 
have 

B„, fc (y + 1)= ^ y |5| f, *J:^ 1 )a„-( t „(S)+i).fc(5n[n-(t n (S) + 2)]) 

SC[n-l] V £n ^ j + 7 
= E E yV 5| " l (^ 1 1 )«„-( l+ i)^(5n[n-(z + 2)]) 



i>0 sc[»-i] 

- i„(S)=i 

Ef' + iV E Mi+iu^nf"-^ 2 )])^ 1 " 

i>0 ^ ' SC|n-l] 

t„(S)=i 

Ef.+JV e «^ +1 u.(^ isi 

j>0 V / SC[n-(i+l)] 

E0+i 1 V B MWk(» + 1 ) 
^f fc + 1 ) y '-%_ i , fe ( y + l). 



i>i 



□ 

Proof of Theorem [71 Replacing n and ?/ by n + + 1 and y — 1, respectively, in 
Proposition Q] yields the recurrence (j4|): 

B n+k +i,k(y) = E ( ■ ) ( y _ 1 ) l_1 ^n+fc+l-i,fe(2/) 
i=l ^ z ' 

for n > 0, with the initial conditions P>i^(y) = Ai(y), < i < k. □ 

Consequently, by multiplying the above recurrence by z n and summing over all 
n > 0, we have the generating function Bfe(z,y): 

1 + £ L(y) - E + ') ^ - l)'- 1 ^^ ** 



i-ECt 1 )^-^ 1 
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3.2. Solving the recurrence for B nJ ±. Before we proceed to solve the recurrence 
for B n: k, we first examine the special case of k = 4 which is quite illuminating. We 
note that the characteristic polynomial for the recurrence for B n ^ is 

h{z) = z 5 - 5z 4 + 10(1 - y)z 3 - 10(1 - y) 2 z 2 + 5(1 - yfz - (1 - y) 4 
= (z - 1 + yf - yz 5 

i-y 

Substituting y — t 5 in the expression above, we see that the roots of h{z) are just 

1 - t 5 

pAt) = -, 0<j<4, 

where uj = cxp(^p) is a primitive 5th root of unity. Hence, the general term for 
B n ,4:(t) can written as 

4 

fl n , 4 (t) = 

t=0 

where the ai(t) are appropriately chosen coefficients (polynomials in t). To deter- 
mine the cti(t) we need to solve the following system of linear equations: 

4 

J2<Xi(t)4(t) = B jA (t) = A At% < j < 4. 

i=0 

Thus, cti(t) can be expressed as the ratio N^i + \(t)/Di(t) of two determinants. The 
denominator £> 4 (i) is just a standard Vandermonde determinant whose + + 
entry is pj(t). The numerator N^i + \(t) is formed from D±{t) by replacing the 
elements p\(t) in the (i + l)st row by Aj(t 5 ). A quick computation (using the 
symbolic computation package Maple) gives: 

D 4 (t) = 25V5(1 -i 5 ) 6 t 10 ; 
JV 4 ,i(t) = 5y/5 (t 12 + t w + 2t 9 + At 8 + At 7 + At 5 + At 4 + 2t 3 + t 2 + 1)(1 - t 5 f(l - tft 10 

and, in general, iV^j+i (i) = N^\(u l t). 

Substituting the value a (t) — N 4 ,i(t) / D '4(f) into the first term in the expansion 
of B ni 4, we get 

a (t)(l + t + t 2 + t 3 +t 4 ) n 
= \{t 12 + 1 10 + 2t 9 + At 8 + 4t 7 + At 5 + At 4 + 2t 3 + 1 2 + 1)(1 + f + 1 2 + t 3 + t 4 ) n - 3 . 

Now, since the other four terms a»(t)(l + i + t 2 + t 3 + t 4 ) n arise by replacing t by 
to l t then in the sum of all five terms, the only powers of t that survive are those 
which have powers which are multiples of 5. Thus, we can conclude that if we write 

(t 12 + i 10 + 2t 9 + At 8 + At 7 + At 5 + At 4 + 2t 3 + 1 2 + 1 ) ( 1 + 1 + 1 2 + 1 3 + 1 4 ) n ~ 3 = ^ /3 (r )t r 

r 

then 6 n , 4 (d) = /3(5d). In other words, the number of permutations tt € B n ^ with d 
descents is given by the coefficient of t 5d in the expansion of the above polynomial. 
Incidentally, the corresponding results for the earlier B n ,i are as follows: b n ,i(d) = 
(3(2d) in the expansion of 

(l + tr=J2/3(r)t r , 

r 

so b nm i(d) = (^); b nm2 {d) — /?(3d) in the expansion of 

(i+t 2 )(i+tr- i = j2mt r ; 
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and b n ^(d) = /3(id) in the expansion of 

(i + i 2 + 2t 3 +t 4 + t 6 )(i + t) n - 2 = J2 /W- 

r 

The preceding arguments have now set the stage for dealing with the general 
case of B Ui k- Of course, the arguments will be somewhat more involved but it is 
hoped that treating the above special case will be a useful guide for the reader. 

3.3. Solving the recurrence for B n ^- 

Theorem 2. We have B„ jk (y) = Y. d f j k{{k + l)d)y {k+1)d , where 

/ 1 fc+i x n ~ k 

Eftc?> J =A(u)(-=^- 



and 

p k (u) = j2 A ^(n k+1 ){u k+1 - iyJ2 

3=0 «=3 



Proof. To solve ^ for we first need to compute the roots of the corresponding 
characteristic polynomial 

z n+k+l _ ^ + ^ (j, - l)-l z »+fe+l-« 

which can be rewritten as 

Substituting y = u k+1 , this becomes 

(z - (1 - u k+1 )) k+1 -u k+1 z k+1 . 
The k + 1 roots of this polynomial are easily seen to be the expressions 

1 - u k+1 

pi = Pi(u) = - — , < i < fc, 

1 — ui l u 

where uj = exp(^pj) is a primitive (fc + l)st root of unity. Thus, we can express 
B ni k in the form 

k 

B n , k (u k+1 ) = $>m(«)p? 

i=0 

for an appropriate choice of coefficients a k ,i (which depend on the initial conditions). 
In fact, writing down the expressions for the first k + 1 B n:k s, we have: 

k 

B 3 ,k{u k+1 ) = X) a *.*(«W = for < j < fc. 

i=0 

We can solve this as a system of k + 1 linear equations in the fc + 1 unknown 
coefficients a kl i(u) by representing the solution in the usual way as a ratio of two 
determinants. In particular, the expression for a k ,o is given by 

deti? fe (u) 
det S k (u) 

where S k (u) and Rk(u) are (fc + 1) by (fc + 1) matrices defined by 
S k (u) = (S k (i + 1,3 + 1)) with S k (i + l,j + 1) = p> 
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and 

' Aj{u k+l ) if i = 0, 
p{ if i > 0. 



Rk(u) = (ifc(t + 1,3 + 1)) with R k (i + l,j + 1) 



Now the bottom determinant is a standard Vandermonde determinant which has 
the value 

detSfc(it)= Yl (ft - / 9 *)- 

0<i<j<k 

The top determinant is almost a Vandermonde determinant (except for the first 
row). Its value has the form 

det R k (u) = Top k (u) ■ (pj - Pi ) 

0<i<j<k 

where Top fc (u) is a polynomial in u which we will soon determine. Hence, in the 
ratio ©, the terms which do not involve po cancel, leaving the reduced form 

/ \ T °Pfc(«) 

otkMu) - 



rij>o(Pj - Pa) 

However, we have 

j>0 j>0 v 

{\-u k+i ) k u k rw 1 -^') 



(i-u) k rw 1 -^) 

(l-u k+1 ) k u k k+1 _{l-u k+1 ) k ~ 1 u k 



(l-u) k {l~u) k 



fc-i 



(k + 1). (7) 



On the other hand, for the top we have by standard properties of Vandermonde 
determinants: 

To Pk (u) = Y l (-l) k - i A k - j (u k+1 )S k , j (p 1 ,..., Pk ), (8) 

j>0 

where S k j(xi, . . . ,x k ) is the elementary symmetric function of degree j in the k 
variables x\,...,x k (and we recall that A t (y) = J2j>o (j)^')- Now consider the 
generating function 

X k (z) = (z- p )(z - Pi) ■ ■ ■ (z - p k ) 

= J2(-^) t S k+ iAPo,---,p k )z k+1 - t 



1 - u 



k +l 



z — 



1 — u % u 



n 

i>0 

U i > (z-(l-n k+1 )-^uz) 

rW 1 -^) 

(z-{l-u k+1 )f +1 -u k+1 z k+1 

1 - ' 
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What we are interested in is 

Y k (z) = = J2(-l)%dPi, • ■ • ,Pk)z k ~ ' 

z - Po ^ 



i>0 



(z - (1 - u k+1 )) k+1 -u k+1 z k+1 



(1- 



1-u 

A- 



— ^(2-l +!1 l+1 )'( M zf- 



1 - U fc +! 

i^ed-i) , (;)( i -« w ) , « m ^ 

j=0 i=j V 



Thus, by identifying the coefficient of z k 3 , we have 



k 



(1 - u)(l - u k+1 y >A /i 

Sfcj(A»i,...,p fc ) = 1 _ Mfc+1 ( 7 - 

Therefore, we have 

Topju) = 2(-l)*-^_ i («^ 1 )S w 0' 1 , . . . ,p fc ) 



J>0 

n ,,fc+i\?" fe /„-\ 



J>0 i=j 



As a consequence, we find by © and ([7]) that 
«fc,o(it)Po 

Top fc ( U ) /l-^V 



(-l) fe n 7 >o(Po-p,) V i- 



(9) 



1 - „k+ 1 \ n ~ k k k r; 

-L U \ X ^ / . s o , k_I_1 >. k_I_1 >. o X ^ / 6 



J=0 i=j 



(fc + l)?i fc V 1-u 

(We should keep in mind that this expression actually is a polynomial in u.) To 
determine the other coefficients a k ,t(u), we make the following observations. First, 
for 1 < t < k, we can cyclically permute the rows of the (Vandermonde) matrix 



S k (u) to form the new matrix Sg\u) = (S { k t} (i + l,j + 1)) = (p{ +t ). We 



can 



then form the corresponding matrix R^\u) by replacing the top row of S^\u) by 
Ao(u k+1 ), . . . , A k (u k+1 ). In this way, we can express the coefficient ctk,t(u) as: 

a 



, . deti^O) , . 

«*.*(«) = . „m ■ ( 10 ) 



detS^(«) 

However, observe that the resulting computations for determining ctk,t(u) are ex- 
actly the same as those for ak,o(u) where we replace u by w t u. This is because 
Aj(u k+1 ) = Aj(u)*u) k+1 . Consequently, 

a k ,t{u) = a kfi (u t u). 
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Therefore, for each n, the expression ctk,t{u)Pt as a polynomial in u can be obtained 
from akfi{u)pQ by replacing u by o/u. However, this implies that in the sum 
2t=n a k,t[ u )pt ■> the only terms that survive are those powers u m of u which are 
multiples of k + 1, since for m ^ (mod fc 4 1), we have 2i=o = ^. One the 
other hand, for u m with m = (mod fc 4 1), we must multiply the coefficients by 
k 4 1 since in this case, all the powers a/ are 1. Consequently, if we write 

(fc + l)a fc ,o(uK = X;&OV (11) 
i 

then we have 

s^(y) = E^(( fc + 1 ) d )y (fe+1)<i - 

In other words, the number 6 ni fc(r) of permutations in B n> k having exactly r descents 
is equal to the fc + 1 times the coefficient of 

If we express © (times fc + 1) in the form 

P k (u) (1 + u 4 ■ ■ ■ 4 u k ) n ~ k (12) 
then the first few values of P^ (u) are shown below. 



fc 


Pk(u) 


1 


1 












2 


1 + U 












3 


1 + m4 


2m 2 H 


- u 3 


+ M 4 






4 


1 +«4 


2u 2 H 


- 4u ; 


' + 4u 4 + 4u 5 - 


F- 4u 6 4 


2u 7 + M 8 + u 9 


5 


1 4u4 


2u 2 H 


- 4u ; 


! + 8u 4 + llu 5 


4 llu 6 


4 14u 7 + 16u 8 4 




+14u 9 


+- llu 


10 + 


llu 11 + 8w 12 H 


- 4u 13 4 


-2m 14 + u 15 + m 16 



There is clearly a lot of structure in the polynomials Pk(u) which will be discussed 
in the next section. 

4. The structure of Pk(u) 
Let us first write down the expression for Pk(u) which came from (9): 

k k / . \ 

^(«)=X;(- l ) , '^-i(«^ 1 )( l -«^ 1 ) i E (!•)""*■ ( 13 ) 

j=0 i=j w 

k 2 

If we write Pk(u) = ^~V_ Q o^it 4 , we will define the stretch of Pfc(tt) to be 

fc — 2 

PP fe ( U ) = a + E E «i +J+ ( fe +i)" 2+J+(fe+lb+J - 
i=o j=a 

What this does to Pk(u) is to insert coefficients at every (fc + l)st term, starting 
after cxq. Thus, the stretched polynomials corresponding to the values of Pk(u) 
given in the array above are: 



fc 


PPk{u) 


1 


1 










2 


1 + u 2 










3 


1 + u 2 4 


2m 3 - 


fll 4 +M 6 






4 


1 + u 2 4 


2m 3 - 


f- 4u 4 + 4u 5 


+ 4m 7 4 


- 4u 8 4 2u 9 4 u 1Q 4 u 12 


5 


1 + u 2 4 


2u 3 - 


f 4u 4 + 8m 5 


+ llu 6 


4 llu 8 4 14u 9 4 16u 10 4 




+ 14u n 


+- 11m 


12 4llu 14 4 


8u 15 4 


4m 16 4 2m 17 4 u 18 4 u 2Q 



Note that if Pjt('u) has degree fc 2 then PPk(u) has degree fc 2 4 fc. 
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Theorem 3. For all k>l, 

P k +i(u) = PP k {u) • (1 + u + u 2 + • • • + u k+1 ). 
Proof. From ([TB")) and the definition of PP k (u), we can write 



fc 

' S 



pp k ( u ) = - u k +yA k ^ t (u k+2 ) J2 

i=0 s=t 

We want to show that 

f - u k+2 

P k+1 (u) = PP k (u) • — . 

1 — u 

Thus, it is enough to prove that A = B where 
A = PP k+1 (u)(l-u k+2 ) 

k k 

= j2(-iy(i-u k+2 r +i A k _ t (u k+2 )j2 

t=0 s—t 

and 

B = P k (u)(l-u) 

fc+l fc+l / \ 

= - u fc+2 )^ fc+1 _ t (u fe + 2 ) x: - u). 



t=0 

Observe that 



fc+l /fc+l /x > 

B = £(-!)*(! - U fe+2 )^ fe+1 - t ( U fe+2 £ U K ^ ~ 



t=0 \s=t 



fe+l /fc+l /X \ 

= $>i) 4 (i ^ u fc+2 )^ fe+1 _ t ( U fe+2 ) fx: - »)) 
+A fe+1 ( M fc+2 )(f-i t -' £ - 2 )(- M ) 

= |> 1)t(1 - ^ +2 ) t+1 ^-^ +2 ) ( E ( t + - ^ 

-fc-2- 



t=0 \s=t+l 
,- k -2 



-u(l-u z )A k+1 

= £(-l)*(l - ^ +2 ) t+1 A,_ t ( M fc + 2 ) ( J] 
=o 

- u(l - u- k - 2 )A k+1 

= 5>i)*(i - u k + 2 y+iA k ^ t (u k + 2 ) [ v , / ) K , + ] 



t=0 

-fc-2\ 



-u{l-u z )A k+1 . 
Hence, to prove that A = B, we only need to establish 

J2(-iY(i-u k + 2 y+iA k _ t (u k + 2 ) (-(^+iy- k -^-u(i- u - k - 2 )A k+i 

= u{l-u- k - 2 )A k+1 {u k+2 ) 

However, this would follow from 

fc' 



53(a:- iyA k ,- t (x)( k ^j = xA k ,{x) (14) 



t=0 
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by taking k! = k + 1 and x = u k+2 . So it remains to prove ((HJ). 

To do this we will use the standard generating function for A n (w) from ((3J): 



Consider 



and 



iffl/ n! ' ™ nl e (w-i)z _ w 

n,m>0 n>0 



fc>0 



fc>0 t=0 

It will suffice to show that F = G. Now 



Ai>0 t>o ^ ' 



££(*-i)*A fc _^)(f)^-i 



t>0 fc>t 



t>0 fc'>0 



vV fc! 
fc' + A z fc ' +t 



i /(fc' + t)! 



t>0 fc'>0 



,(x-l), . jjlf - 1 = < l " 



3 (a;— 1)^ <£ ^ l)z 



On the other hand, 



F(x,z) = ^xA k (x 



fc! 

fc>0 



l-x \ _ x(l -e^- x >) 



g(x — l)z j. J 1)5 

as desired. This completes the proof of Theorem [3] □ 
Theorem 4. The coefficients of P k (u) are symmetric and unimodal. 

Proof. It follows from Theorem [3] that we can construct the coefficient sequence 
for Pfe_|_r(tt) from that of Pk{u) by the following rule (where we assume that all 
coefficients of u* in Pk{u) are if t < or t > fc 2 ). Namely, suppose we write 

Pk(u) — 2i=o a i ul so that we have the coefficient sequence A k — (oto, oti, ■ ■ ■ , ct^). 
Now form the new sequence B k = (/?o, Pi, ■ ■ ■ Pk 2 +k) by the rule 

i 

ft = E a J j - * - fc2 + k - 

j=i-k 

Finally, starting with /3q, insert duplicate values for the coefficients 

/3o)/3fc+i)^2(/b+i)) • • • )A(fc+i)j • • • )/3(fc-i)(fc+i) and /3fe(fe + i). 
Thus, this will generate the sequence 

(Po, Po, Pi, P2, ■ ■ ■ , Pk, Pk+i, Pk+i, Pk+2, ■■ ■ , Pk 2 +k-l, Pk 2 +k, Pk 2 +k)- 



DESCENT POLYNOMIALS FOR PERMUTATIONS WITH BOUNDED DROP 15 

This new sequence will in fact just be the coefficient sequence A^+i for Pk+\{u). 
For example, starting with Pi(u) = l + u,we have A\ = (1, 1) and so B\ = (1, 2, 1). 
Now, inserting the duplicate values for /3q = 1 and fa = 1, we get the coefficient 
sequence A 2 = (1, 1, 2, 1, 1) for P2{u) = 1 + u + 2u 2 + u 3 + u 4 . Repeating this 
process for A 2 , we sum blocks of length 3 to get B 2 = (1,2,4,4,4,2, 1). Inserting 
duplicates for entries at positions 0, 3 and 6 gives us the new coefficient sequence 
A 3 = (1, 1, 2, 4, 4, 4, 4, 2, 1, 1) of P 3 = l+u+2u 2 +4u 3 +4u 4 +4u 5 +4u 6 +2u 7 +u*+u 9 , 
etc. It is also clear from this procedure that if Ak is symmetric and unimodal, then 
so is -Bfe, and consequently, so is A^+i- This is what we claimed. □ 

4.1. An Eulerian identity. Note that since Pk(u) is symmetric and has degree 
u k , we have Pk(u) — u k Pfc(~)- Replacing Pk(u) by its expression in (fll?|) . we 
obtain (with some calculation) the interesting identity 

jy-iy Q (i- x yA a+b ^( x ) = X J2 Q (i-xyA a+b ^(x)+( a b + 6 ) (i- x y+ b+1 

for all integers a and b provided that a + b > 0. 

References 

J. Buhler, D. Eisenberg, R. Graham and C. Wright, Juggling drops and descents, Amer. 
Math. Monthly 101 (1994), 507-519. 

F. Chung, and R. L. Graham, Primitive juggling sequences, Amer. Math. Monthly 115, 
March 2008, 185-194. 

L. Euler, Methodus universalis series summandi ulterius promota, Commentarii academiae 
scientiarum imperialis Petropolitanae 8 (1736), 147-158. Reprinted in his Pera Omnia, series 
1, volume 14, 124-137. 

D. Foata and G.-N. Han, q-series in Combinatorics; permutation statistics (Lecture Notes), 
preliminary edition, 2004. 

D. Foata and G.-N. Han, Fix-Mahonian calculus III: A quadruple distribution, Monatsh. 
Math. 154 (2008), no. 3, 177-197. 

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison- Wesley, 1994. 
D. E. Knuth: The Art of Computer Programming, Vol. 1, Fundamental algorithms. Addison- 
Wesley, Reading, 1969. 

D. E. Knuth, The Art of Computer Programming, Vol. 3, Sorting and Searching, Addison- 
Wesley, Reading, 2nd cd. 1998. 

P. A. MacMahon, Combinatory Analysis, 2 volumes, Cambridge University Press, London, 
1915-1916. Reprinted by Chelsea, New York, 1960. 

O. Rodrigues, Note sur les inversions, ou derangements produits dans les permutations, J. 
de Math. 4 (1839), 236-240. 

J. Sharcshian and M. L. Wachs, g-Eulerian polynomials: exccdance number and major index. 
Electron. Res. Announc. Amer. Math. Soc. 13 (2007), 33-45. 

J. Shareshian and M. L. Wachs, Eulerian quasisymmetric functions, preprint 2009. 



[1 
[2 
[3 

[1 
[5 

[e; 
\- 

[9 

[io; 
in 

[12 



F. Chung and R. Graham: University of California at San Diego, La Jolla CA 92093, 
USA 

A. Claesson: The Mathematics Institute, School of Computer Science, Reykjavik 
University, 103 Reykjavik, Iceland 

M. Dukes: Science Institute, University of Iceland, 107 Reykjavik, Iceland 



